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Abstract 

It is long known that the Fokker-Planck equation with prescribed constant 
coefficients of diffusion and linear friction describes the ensemble average of 
the stochastic evolutions in velocity space of a Brownian test particle immersed 
in a heat bath of fixed temperature. Apparently, it is not so well known that 
the same partial differential equation, but now with constant coefficients which 
are functionals of the solution itself rather than being prescribed, describes 
the kinetic evolution (in the A —> oo limit) of an isolated A-particle system 
with certain stochastic interactions. Here we discuss in detail this recently 
discovered interpretation. 

KEYWORDS: Kinetic theory, Kac program, propagation of chaos, diffusion equa¬ 
tion on a high-dimensional sphere, Fokker-Planck equation. 


1 Introduction 


As is well known, luhOrfinl . imi . , the ensemble average of the stochastic 

evolutions in velocity space of a Brownian test particle^ of unit mass, immersed in 
a drifting uniform heat bath of fixed temperature T and constant drift velocity m, 
is governed by the Fokker-Planck equation with prescribed constant coefficients of 
diffusion and (linear) friction, 

dtf{v\t) = ■ (rd^fiv-t) + {v- u)f{v,t)y (1) 

Here, /(. ;t) : ^ K+ is the ensemble’s probability density function on velocity 

space at time t S ]R_|_, and an overall constant has been absorbed in the time 
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^For the beginnings of the theory of Brownian motion, see the collection of Einstein’s papers 
with commentary IEin85l . 


1 










scale. Of course, we could also shift v to obtain u = 0, then rescale v, t, and 
/ to obtain T = 1; however, for pedagogical purposes we refrain from doing so. 
The solution of is given by f{v;t) = Gt{w,v\u;T)fo{w)d^w, where 

foiv) = f{v;0) and 


Gt(u;, u|m; T) = (27rT(l — e ^‘)) ^ exp 


/ 1 \v — u — we \ 

V ^ 1 - e-2* J 


( 2 ) 


is the Green function for o , see luhOr.ml . ITMI . In its standard form, i.e. with 
r = 1 and M = 0, lI21) is known as the (Mehler) kernel of the adjoint Ornstein- 
Uhlenbeck semigroup (a.k.a. Fokker-Planck semigroup). 

Over the years, the Ornstein-Uhlenbeck semigroup and its adjoint have come to 
play an important role in several branches of probability theory |HSU99| related, 
in some form, to Brownian motions. The fact that the explicitly known kernel II 
of the Fokker-Planck semigroup readily lends itself to analytical estimates has led 
to useful applications also outside the realm of probability theory. In particular, in 
recent years the Fokker-Planck semigroup has found applications in kinetic theory, 
the subheld of transport theory which is concerned with the approach to equilibrium 
and the response to driving external forces of individual continuum systems not in 
local thermal equilibrium; see, for instance, the review 

However, the linear Fokker-Planck equation itself, m, usually is not thought 
of as a kinetic equation for the particle density function on velocity space of an in¬ 
dividual, isolated space-homogeneous system of particles in some compact domain, 
which perform a microscopic autonomous dynamics that may be deterministic or 
stochastic but should satisfy the usual conservation laws of mass (particle num¬ 
ber), energy and, depending on the shape of the domain in physical space and its 
boundary conditions, also momentum and angular momentum. Evidently the very 
meaning of / and the parameters u and T in o voids this interpretation. Yet, with 
a re-interpretation of /, u and T it is possible to assign to o a kinetic meaning. 

Incidentally, the hrst result showing that at least a partial re-interpretation of 
0 in this direction is possible can be found in a paper by Villani mm who. 
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in his study of the space-homogeneous Landau equation for the weak deflection 
(i.e. Landau) limit of a gas of particles with Maxwellian molecular interactions, 
discovered that for isotropic velocity distribution functions / (and only for these) 
the Landau equation is identical to o, with parameters u = 0 and T matched to 
guarantee energy conservation. For general non-isotropic data the Landau equation 
for Maxwell molecules is identical to a more complicated equation than O- 

To pave the ground for a complete re-interpretation of o, which requires re¬ 
assigning the meaning of /, u and T, we first note that by the linearity of o we 
can scale / to any positive normalization we want. We now introduce the following 
functionals of /, 
the “mass of /” 

m{f)= [ f{v\t)d^v, (3) 

the “momentum of /” 

P{f)= [ vf{v;t)d^v, (4) 

dR3 

and the “energy of /” 

e(/)= / hv\^fiv]t)d^v . (5) 

Js.3 l 

The “angular momentum of /” for a space-homogeneous f{v;t) is simply j{f) = 
a^cM X p(/), with xcM the center of mass of the system, but this does not add any 
further insight and hence will not be considered explicitly. The functionals 0, 0, 
and 0 inherit some time dependence from the solution /(. ;t) of 0, but to find 
this dependence explicitly it is not necessary to solve for / first. Indeed, it is an 
elementary exercise in integration by parts to extract from o the following linear 
evolution equations with constant coefficients for m, p, and e, 

m = 0 , (6) 

p = mu — p, (7) 


3 


e = 3T — 2e + M • p, 


( 8 ) 
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which, beside the conservation of mass, i.e. m{f) = to(/o), describe the ex¬ 
ponentially fast convergence to a stationary state p{f) m{fo)u and e(/) 

|T -|- im(/o)|Mp. While all this is of course quite trivial and well known, the rel¬ 
evant fact to realize here is that whenever the energy and the momentum of the 
initial fo equal these asymptotically stationary values, viz. if p(/o) = m(fo)u and 
e(/o) = + ^to(/o)|mP, then beside mass m, also energy e and momentum p will 

be conserved. Conservation of mass, energy, and momentum for such a large subset 
of initial data fo does not yet mean that we may already think of the linear equa¬ 
tion o as a kinetic equation, which should conserve mass, energy, and (depending 
on the shape of the domain in physical space and its boundary conditions) also 
momentum for all initial data, no matter what their mass, energy and momentum 
are; moreover, a genuine kinetic equation for particles with (pair or higher order) 
interactions must express the time derivative of / in terms of an at least^ bilinear 
operator in /. However, with the help of ®, 0 and we now replace T and u 
in o to obtain just such a kinetic equation. 

Indeed, consider the a priori nonlinear Fokker-Planck equation 

dtf{v] t)=d^- (2e(/)m(/) - |p(/)P)9„/(u; t) + {m{f)v-p{f))f{v, t)^, (9) 

where /(.; t) : K.^ —> R+ now is a particle density function on velocity space at time 
t G M+. The right-hand side of 0 is a sum of a bilinear and a trilinear operator 
acting on / which now guarantees conservation of mass, momentum, and energy 
for all initial data fo ^ 0, as verihed by repeating the easy exercise in elementary 
integrations by parts using @ to find m = 0 as well as p = mp — pm = 0 and 
e = 2em— \p\'^m — 2em+ |ppm = 0. Of course, after this fact of mass, momentum, 
and energy conservations the a priori nonlinear equation © in effect becomes just 
a completely and explicitly solvable linear^ Fokker-Planck equation 0, only now 
with parameters u and T which are not prescribed but determined through the 

^The Boltzmann, the Landau, and the Vlasov kinetic equations have bilinear “interaction 
operators,” the Balescu-Lenard—Guernsey equation has a higher order nonlinearity which reduces 
to the bilinear format in the long wavelength regime. 

^In this sense ini is “almost nonlinear,” or “essentially linear,” depending on one’s viewpoint. 
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initial data fo, viz. u = p{fo)/m{fo) = Uq and |T = e(/o) - |p(/o)P/2m(/o) = Eq; 
we also set m(/o) = mg and e(/o) = eg. Accordingly, (0 inherits from the 
feature that, as t —> oo^ its solutions / converge pointwise exponentially fast to the 
Maxwellian equilibrium state 



( 10 ) 


with monotonically increasing relative entropy 



( 11 ) 


which in fact approaches its maximum value 0 exponentially fast. 

Since 0 displays all the familiar features of a kinetic equation (formal nonlin¬ 


earity; conservation laws of mass, momentum, energy; an i/-Theorem; approach to 
equilibrium; Maxwellian equilibrium states), at this point we may legitimately con¬ 
template as a kinetic equation of some spatially homogeneous, isolated system 
of N interacting particles in a compact spatial domain compatible with momentum 
conservation (e.g. a rectangle with periodic boundary conditions). In the remain¬ 
der of this paper we show explicitly how m arises from the Kolmogorov equation^ 
for the adjoint evolution of an underlying A^-particle Markov process in the limit 
N ^ oo. We use the strategy originally introduced by Kac |Kac5fi| in 1956 in 
the context of his work on a caricature of the Boltzmann equation; for important 
recent work on Kac’s original program, see 1^^1.921 . As Kac realized, the crucial 
property that needs to be established in order to validate the N oo limit is what 
he called “propagation of chaos,” which loosely speaking means that if the particle 
velocities are uncorrelated at t = 0, they remain uncorrelated at later times; this 
can be rigorously true only on the continuum scale in the limit N —>■ oo. 

Interestingly enough, by adding some suitable lower order terms to the puta¬ 
tively simplest fV-particle Markov process that leads to the (kinetic) Fokker-Planck 
equation in the limit N —> oo, the corresponding Kolmogorov equation for an en¬ 
semble of such isolated iV-particle systems can be simplified to be just the diffusion 

^In the physics literature, the Kolmogorov equation for an Af-particle Markov process is tradi¬ 
tionally called “master equation”. 
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equation on the 3-/V — 4-dimensional manifold (a sphere) of constant energy and mo¬ 
mentum. Since therefore both the finite-and the infinite-iV equations are exactly 
solvable, the kinetic limit N oo can be carried out explicitly and studied in great 
detail. For this reason we actually defer the discussion of the underlying 7V-particle 
process to Appendix Ab while in the main part of our paper we analyse the diffusion 
equation on and derive from it the kinetic Fokker-Planck equation on 

Technically, we apply the Laplace-Beltrami operator to a probability density on 
and then integrate out N — n velocities over their constrained domain of 
accessibility. Taking next the limit N ^ oo yields a Fokker-Planck operator act¬ 
ing on the n-th marginal density on Thus we obtain a linear Fokker-Planck 

hierarchy of equations indexed by n. Using the Hewitt-Savage decomposition the¬ 
orem, the hierarchy is seen to be generated by the single, a priori nonlinear kinetic 
Fokker-Planck equation m which in view of the conservation laws is equivalent to 
the essentially linear Fokker-Planck equation o with constant parameters which 
are determined by the initial data. 

Experts in probability theory may have noticed a similarity between the first 
part of our program and what has been called the “Poincare limit” |Bak02| : in fact, 
our approach is “dual” to Bakry’s approach. More specifically, Bakry |Bak02| has 
shown that the action of the Laplace-Beltrami operator for ^ R'^+^ on a 

probability density function over a “radial” coordinate axis of becomes iden¬ 
tical, in the limit N —^oo,to the action of the Ornstein-Uhlenbeck operator on 
the same density viewed as a function over R. Obviously, whenever the “radial” 
function is obtained by taking the marginal of a probability density over i-C. 

by integrating out the — 1 Cartesian coordinates of the embedding space which 
are perpendicular to a fixed “radial” direction, the Ornstein-Uhlenbeck operator 
acts on the limiting marginal density as N oo. This relationship between the 
operators is reflected at the spectral level by the convergence of the whole structure 
of orthogonal eigenfunctions of the Laplacian on (hyper-spherical harmonics) 
to the orthogonal eigenfunctions of the Ornstein-Uhlenbeck operator on R (Hermite 
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polynomials multiplied by the square root of their Gaussian weight function); one 
of the earliest works is | Meh66| . while more recent works on the Poincare limit, con¬ 
taining interesting connections with the theory of Markov semigroups, are |Bak02| 
and |BaMa,n8] . Our procedure is “dual” to Bakry’s approach in the sense that 
we integrate out subsets of the Cartesian variables of the embedding space after 
having applied the Laplace-Beltrami operator to a probability density on the high¬ 
dimensional sphere, thereby obtaining the adjoint Ornstein-Uhlenbeck operator 
acting on the respective marginals; in addition, while Bakry considers only mass 
and energy conservation, we consider conservation of mass, energy, and momentum. 

Incidentally, our work is not inspired by Bakry’s works on the Poincare limit, 
nor by Villani’s discovery about the isotropic evolution of the space-homogeneous 
Landau equation, about both of which we learned only after our own findings. 
Rather, our study of the diffusion equations on the 3N — C-dimensional spheres 
of constant energy (C = 1), respectively energy and momentum (G = 4), which 
began in |KiLa04 j . was originally conceived of as a technically simpler primer for 
our investigation (also in |KiLa04| l of the Balescu-Prigogine master equation for 
Landau’s kinetic equation. And while the present paper is also a technical contin¬ 
uation of |KiLa04| . in the sense that here we supply various calculations that we 
had announced in |KiLa04 |. the main purpose of the present paper is to amplify the 
conceptual spin-off of our technical investigations, the new physical interpretation 
of one of the simplest and best known linear transport equations as an (almost 
nonlinear) kinetic equation. As should be clear from our discussion in this intro¬ 
duction, this kinetic theory interpretation of the prototype Fokker-Planck equation 
may have been suspected by others long ago, yet we have not been able to find the 
whole story in the literature. 

In what follows, for the sake of simplicity we set toq = 1, and accordingly® obtain 
p(/o) = Ug and e(/o) — |p(/o)P/2 = eg — \ug\^/2 = Eq- With these simplifications 

® Setting mo = 1 means we should now speak of the energy per particle eo, the thermal energy 
per particle eo, S'Ud the momentum per particle Pq{= ug). 

















0 now becomes 


dtf{v; i) = • (^^Sod^f{v-1) + (v - Uo)fiv, t)^ . (12) 

While 11211 is essentially a linear PDE, it should just be kept in mind that eo and 
Uq are functionals of / which are determined by the initial data /o and not chosen 
independently.® We next shall derive m from the diffusion equation equation on 
the spirit of Kac’s program. 

2 The Finite-Ensembles 

Consider an infinite ensemble of i.i.d. random vectors where each V = 

(vi, ...jVjv) S represents a possible micro-state of an individual system of N 
particles with velocities Vi = {vn^va^va) G K® and particle positions assumed to 
be uniformly distributed over a periodic box; hence, particle positions will not be 
considered explicitly. Each V takes values in the 3N — 4-dimensional manifold of 
constant energy eo and momentum Uq, 

N N 

= Nuo, - \vkf = Neo, eo > (13) 

k^l k^l 

The manifold is identical to a 3iV — 4-dimensional sphere of radius ^/2N£q 

(where eo appears above llT^ l. centered dX U = {uq, ...,Uo) and embedded in the 
3(JV — l)-dimensional affine linear subspace of given hy U -\- where 

= ]^3Af p g . 'Yk=\'^k = O} is the space of velocities in any 

center-of-mass frame. The ensemble at time r is characterized by a probability 
density (V; r) on the evolution of which is determined by the diffusion 

equation 

a.FW(y;r) = A„3.-4FW(y;r), (14) 

CAQ ,eo 

where A„3Ar-4 is the Laplace-Beltrami operator on Since all particles are 

l'Wl-u,Q,eQ 

of the same kind, we consider only solutions to da which are invariant under the 

®The identification of e with is valid only for isolated systems that can freely translate. 
If a driving external force field F is applied, then e(/) and p(/) are no longer constant and lO — 
with the addition of the forcing term — F ■ dvf its r.h.s. — is the relevant equation. 
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symmetric group Sn applied to the N components in oiV. Clearly, permutation 

symmetry is preserved by the evolution/ We will show that the diffusion equation 
d, here viewed as a master equation, leads precisely to the essentially linear 
Fokker-Planck equation d in the sense of Kac’s program: (a) the Fokker-Planck 
equation m arises as the N ^ oo limit of the equation for the first marginal of 
derived from (ICTl . and (b) propagation of chaos holds. In this section 
we prepare the ground by discussing the finite-fV equation d- The limit N oo 
is carried out in the next section, while propagation of chaos is discussed in the final 
section. 

For the sake of completeness, we begin by listing some general facts about the 
diffusion equation. We note that the Laplacian A„3w-4 is a positive semi-definite, 
essentially self-adjoint operator on the dense domain C 

thus it has a unique self-adjoint extension with domain Its self- 

adjoint extension is the generator of a non-expansive semigroup on 
which is strictly contracting on the orthogonal complement of the constant func¬ 
tions. Thus, we may ask that the initial condition limtio ;t) = Fq^\.) G 

(which implies Fq^'* G Yet, as is well-known, the dif¬ 

fusion semigroup is so strongly regularizing that we may even take (.) G 
probability measure, and obtain F(^)(.;r) G for 

all r > 0. 

In fact, the solutions of d can be computed quite explicitly in terms of an 
eigenfunction expansion. Since via translation by U (choosing a center-of-mass 
frame) and scaling by \/2Neo (choosing a convenient unit of energy) the manifold 
can be identified with the unit sphere centered at the origin of the linear 
subspace C the complete spectrum of A lyrSAi—4 Rnd ctn ortli 020 ii&l bcisis 

iV^-uO,eo 

of eigenfunctions can be obtained from the well-known eigenvalues and eigenfunc- 
tions for the Laplacian on the unit sphere §37V-4 ^ r 3^“3. Of course, in our 

^In what follows, for the sake of notational simplicity we will not enforce this symmetry explic¬ 
itly, but the reader should be aware that (for instance) all the eigenfunctions that appear below 
in the solution for can be easily symmetrized. 
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case the embedding is ^ with isomorphic by a rotation to 

standard Thus we start from and we first carry out a rotation in 

that transforms V to W = UV in such a way that is mapped to the 

ZN — 3-dimensional linear subspace {W : wn = O}- Obviously, is the linear 

transformation that diagonalizes the projection operator onto A complete 

orthonormal set of eigenvectors for such a projection is readily calculated and leads 
to 


Wi 


Wn 


/TV- 1 

Vl 

TV 



N — n 


N 


-wVtE»‘ 


i=2 


N — n + \ 


1 

N — n 


N 




WN-1 = 


Wn = 


V2 


[uat-I - VAr] 


1 ^ 


(15) 


It is easily checked that the matrix associated with this transformation is indeed 
orthogonal, and that wn vanishes whenever V G More generally, the affine 

subspace U is mapped to the linear manifold {W : wn = VNuo} and 

is mapped to 


IW : WN = VNuo, Y. 1^*1" = ^ l“ol" = \ (16) 

which implies that the truncated vector {wi,... ,wn-i) belongs to the sphere 
^ (in iCfc-coordinates). Thus, the transform U allows one to anal¬ 

yse the A^-particle system with energy and momentum conservation (“periodic box” 
setup) in terms of an (TV — 1 (-particle system with only energy conservation (a “con¬ 
tainer with reflecting walls” setup).® For future reference, we also observe that for 
n fixed and N ^ oo the effect of U reduces to a translation of each of the n ve- 

®The gas in such a container was discussed in our earlier work IKiLa04l . but without detailed 
calculations. Our calculations with the w variables here now supply the relevant details. 
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locities by Uq, in the following sense. Consider a consistent hierarchy of vectors 
of increasing size N, in which lower-iV vectors can be obtained from the higher-N 
ones by truncation (i.e. projection). Suppose that the vectors belong to U 
for all N, apply the transformation in US) and look at the n-th component. Since 
J2^=n+i = -^“0 - Sfci where is independent of N, we find 


lim Wn=Vn- Uq. (17) 

N^OC 


We now recall that the Laplacian is invariant under Euclidean transformations. 
Thus, under our orthogonal transformation 7/, the Laplacian A„3iv-4 becomes the 
Laplacian on in space of truncated vectors ..., u;Ar_i) 

(which will also be denoted by W, at the price of abusing the notation). Since 
Ag3Ar-4 = and the Laplacian on the unit sphere §3Af-4 spectrum 




j(j + 3A — 5), j = 0,1,..., the spectrum of Aj^3iv-4 is 
dj) _j(j+3A-5) 


A 


T3iV-4 


2Neo 


j = 0,1,... . 


(18) 


The eigenspace on ^ for the j-th eigenvalue has dimension 


A(j,3A-3) 


(3A-5-b2j)(3A-6-f j)! 
j!(3A-5)! 


(19) 


and is spanned by an orthogonal basis of hyper-spherical harmonics® on g3JV-4 ^ 
]^3Af-3 order j, here denoted 3A—3), with ^ G Dj = {1,... ,Af{j, 3A —3)} 

and with uj G §3Ar-4^ rpj^g indexing of our Yj^g{Lj;3N — 3) follows the convention 
of |Mul98| for his Yjy and differs from what might have been anticipated from the 
familiar convention for spherical harmonics on S^. Our reason for using tildes atop 
the function symbols is to remind the reader that we will use a normalization of 
the Yj^i{u};3N — 3) which conveniently suits our purposes and does not seem to 
agree with any of the existing conventions, such as in |Mul98| or for the spherical 
harmonics on S^. Our convention is motivated by the analysis of the large N 
behavior of the eigenfunctions, carried out in Appendix B. 

®The hyper-spherical harmonics on S'* are restrictions to S’* C M’*+^ of homogeneous harmonic 
polynomials in R"+^. For j > 0 the restriction has to be non-constant, since lo,l = const.. 
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Hence, the eigenspace of Aj^3n- 4 associated with the j-th eigenvalue in (EJ is 
spanned by the eigenfunctions (W/^/2Neo; SiV — 3), £ G Dj, where W is given 

by EJ for n = 1,..., N — 1. To shorten the notation we introduce 

(V) = |m3^-4|-i ~ ^ ; 3iV - 3 ) ; (20) 


here, the factor is introduced for later convenience. 

In terms of the eigenfunctions G^^^(V), the solution to equation itTHl is simply 
given by the generalized Fourier series 

JU+3N-5) 




-4 

'UQ,eQ 


-1 


jeti i&Oj 


2Neo 


( 21 ) 


with Fourier coefhcients given by 


p(N) _ 

^ 3 / - 




( 22 ) 


where (.|.) denotes the inner product in -C^(M^^gg^). Notice, though, that the 
numerator |g|^^) can be extended to mean the canonical pairing of the G^-'^^s 


with an element of their dual space, which allows us to take F^^'^ to be a measure. In 
particular, we may take to be the Dirac measure concentrated at any particular 
point of The formula 12111 then describes the fundamental solution of 

the diffusion equation da- In any event, whatever Fq^\ makes it evident 
that when r —> 00 the ensemble probability density function on decays 


-41 

1*0,eo I 


33iV-4 




exponentially fast to the uniform probability density 
Fqj'^Gq^i (V), which is the constant eigenfunction corresponding to the smallest 
non-degenerate eigenvalue 0 of the Laplacian. 


3 Evolution of the Marginals 

To study the limit N —>■ 00 for the time-evolution of the ensemble measure, we need 
to consider the hierarchy of n-velocity marginal distributions 

FH^)(«i,...,«„;r)= [ FW(y;T)du„+i...dujv (23) 

/o3(N-n)-4 

»‘u0 ,eo 
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where riuo^o^ is given by all the (vn+i, ■ ■ ■ ,vn) such that 

N n N n 

t)fc = nuq - |■^;fe|^ = 2 iveo-y^|^’fc|^ (24) 

i—1 2 =n+l i—1 

and has domain {(vi,..., t)„) : X]fe=i ^ 4(A^ — n)£o} C 

The evolution equation for the n-th marginal {vi, ..., t>„; r) is obtained by 

integrating m over (v„+i,..., i)Ar) G using the representation of the 

Laplace-Beltrami operator given in (gnii of Appendix Aa. Then, a straightforward 
calculation (previously presented in |KiLa04| l shows that satisfies 

o mn) ^ V a a:a_1v V 


1 




A ij^\ ^'^^kdVjk 


d 


2Ne, H as; ■ (<”• - “»> - "») 

®j'=l 


gp{n\N) 


dvj 


y O f u^yn\N)\_ 

2enN Z^ dv, ^ j 


_d 

dv^ 


(25) 


2£oiV 

Clearly, to obtain the solutions of these equations it is advisable to integrate 
the series solution for F(^)(V;r), |2U. For this purpose, it will be convenient 
to calculate the marginals in terms of the rotated variables W. Changing the 
integration variables^*^ gives 


F(”l^)(t;i ,... ,Vn]T) = ^ At j F^^\V\T)dWn+i ■ ■-dwN-i (26) 

where the integral is over §3(JV^-n)-4- abused the notation {V ; r) 

•y'2Afeo-Z;?=il’"li 

by applying it to what is now regarded as a function of (vi, d„, Wn+i, 

To obtain the series solution for F*^"I^^(V;t) we need to express m in the vari¬ 
ables {vi, ...,Vn,Wn+i, ...,wn-i) and then integrate term by term in the spirit of 
(P|l. To accomplish this we need to choose explicitly a basis of spherical harmonics 
Yj^i on §3^-4^ convenient to do this in an iterative fashion, by assuming that 
a basis is known for the spherical harmonics with one independent variable less, 
here Yk^m{^ 3 N- 5 ',^N — 4) with UJ 3 N -5 G Then, the desired basis is ob¬ 

tained |Miil98| by taking all the elements in the given lower-dimensional basis and 

^^Note that defines a one-to-one linear map with determinant ^ between 

.. .ii]v) and (m„+i,...m]v-l.2v), where zj^ = wj^ - 
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multiplying them by associated Legendre functions of the “extra” variable. In our 
case the (3iV — 3)-th variable will be Wn/^j 2 Neo, the first component of W/^j2Neo, 
and will be a unit vector in the space of the remaining ZN — 4 components, 

denoted by {W)zN-i/\/‘^Nea\ thus, 


K,. 


/ W 

V ’ 


37V - 3 = Ffc 


(^)3jV-4 . 

\/2^eo 


3A^ - 4 P: 






37V - 3 (27) 


where k = 0,1,... ,j, m = 1,... ,JV{k, 37V — 4) and each choice of the pair k, m is 
associated with a value of the degeneracy index (. for the basis Yj^f, moreover, 
is an associated Legendre function |Mii198| . suitably normalized (see Appendix B). 
By repeating this process 3n times, we write out the eigenfunctions in the form 




W 


:;37V-3 =yfe3„, 


{W)3N —3n—3 
\/‘2.Neo 


37V - 3n - 3 x 


(28) 


,pk3n 


( 


Wn3 


fe3n-i\^^2Areo’ 


37V - 3n - 2 P, 


-)k2 


(- 


Wl2 


-A 


■P‘ 


Wii 


_ 

V 


;37V- 3 


V y/2Neo ’ 

where 0 ^ fc 3 „ ^ < fci ^ j and m = 1,... ,N{k 3 n, 37V — 3n — 3). Now let 

... ,Vn) denote the n-th “marginal” of Gj^\v) (as for in (EHll l. 

and set TV* = N — n — 1. We find 




ff3W-4| 


^3 


f jW) 


37V* 


V V2Areo 


;37V* dWn+1 . ■ - dwN-i 


N 

N-n^ fc; 


pksn 


Wn3 


;37V-3n-2 


Wll 


( - _ 

^ V 3/2Neo 


;37V-3 (29) 


^\y^2Neo 

where the integral is over the same domain as in (US). The integral of is non¬ 

zero if and only if fc 3 „ = 0 and m = 1, and the integrals over Yq,! are determined 
only up to the overall factor Yq.Ij which we may choose to be unity without loss 
of generality. Accordingly, = 0 unless £ G "Dj C Dj, where Dj contains 

the indices of the basis functions that “descend” from the uniform harmonic in 
f-3n-3_ Pqj. gy(;.p I’g |;Pg integrated eigenfunctions then become 




Wll 


;3iv-3)...p°3_, 


Wn3 


\,y2Neo' 


3N-3n-2 


§3(7V-n)-4| 

1 

1 

(a 

\ N-n |S3^-4| 

! -3n 

sj2Neo 



3(N-n)-4 


( 30 ) 
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The series for the n-th marginal (■ i "t) (the integrated eu) is a series in the 

functions , viz. 




Vn\T) 


jGNU{0} fgUj 


j{j+3N-5) _ 

Vn)e (31) 


4 The Limit N ^ oo 


We are now ready to take the infinitely many particles limit. First of all, we observe 
that the evolution equation for the marginal velocity densities (i;i,..., Vn, t) = 
limAT^oo ,Vn',T) which obtains in the formal limit N —>■ oo from it^ 

is the essentially linear Fokker-Planck equation in 




(32) 


dv^ ' 2so' 

i — ± 

We now show that the series expansion for the time-evolved finite-iV marginals 
pin\N)^ ^ converge under natural conditions to solutions of these equations. 
Beginning with the spectrum of Ajyj3w-4, we note that the limit N oo yields 


lim|AW r =ll^l 

AT^oo I M«o.«o j j=0 I j 


'j=0 


(33) 


Thus, the limit spectrum is discrete. In particular, there is a spectral gap separating 
the origin from the rest of the spectrum. As a result, the time evolution of the limit 
N oo continues to approach a stationary state exponentially fast when t oo. 
Coming to the eigenfunctions, the expression on the second line in contains 


1-1 


(the j = 0 


the n-velocity marginal distribution of the uniform density |M^^“ 
case). As is well-known at least since the time of Boltzmann, this distribution 
converges pointwise when A ^ oo to the n-velocity drifting Maxwellian on R^”, 



(recall (ITTIl h In terms of eigenfunctions this means that the “projection” onto R^" 
of the j = 0 eigenfunction of the Laplace-Beltrami operator on converges 

pointwise (in fact, even uniformly) to the j = 0 eigenfunction of the linear Fokker- 
Planck operator in R^", appearing in the r.h.s. of The connection between 
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the eigenfunctions generalizes to the cases j 0; cf. |BaMa03] for the special 
case Mq = 0. The asymptotic behavior for ^ oo of the associated Legendre 
functions in which is discussed in Appendix B, together with yields that 
, Vn) = limAT^oo , «„) exists pointwise for all (ui,..., Vn) G 

with 



for all £ G Dj, where Hm{x) is the Hermite polynomial of degree m on R, and 
we defined mi = j — fci,m 2 = ki — k 2 , ■ ■ ■ = fcsn-i- In terms of the mi’s, 

the index set counts all the choices of integers 0 ^ mi,, m 3 „ ^ j such that 
Si=i = j- For n = 1 one readily recognizes the well-known eigenfunctions 
| IIis96| for the linear Fokker-Planck operator in R^, viz. r.h.s.lO with constant eo 
and Mq, easily calculated by separation of variables. In fact, what we have recovered 
are precisely the eigenfunctions for the linear Fokker-Planck operator in R^", see 

(| 23 ). 

Now assume that one can choose sequences of initial conditions such that, 
for each fixed j and I, the Fourier coefficients F^^'^ converge to a limit Fj^f, such 
that each initial n-velocity marginal density, n G N, converges in (£^ n £^)(R^"') to 

/(’")(ui,...,u„;0) = /®”(vi,...,t;„) + ^ (36) 

it then follows that the subsequent evolution of the n-velocity marginal densities is 
given by 

.. ,Vn;T) = f^^{vi,...,Vn)+Y Y (■»!. • • •: 'Wn)e" . (37) 

Formula describes an exponentially fast approach to equilibrium in the en¬ 
semble of infinite systems. The /*-”^(.; t) G (£^ n £^)(R^"), and in addition they 
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automatically satisfy 



f f'^'^\vi,...,Vn,T)dVi...dVr, = l 

jRSn 

/ (vi H- \- Vn)f^""\vi, . . . ,V„]t) dvi... dVn=nuo 

jR3n 

H - h ...,Vn;T)dvi... dVn=neo 


(38) 

(39) 

(40) 


for all T > 0 (recall that cq = Eq+I^^Io/^)- In fact, (HTTIl solves which now implies 
that /^"^(. ;t) can also be expressed through integration of the initial data against 
the n-fold tensor product of (|21). The upshot is that ■ ;t) G S(R^") Vt > 0 
(Schwartz space). To vindicate these conclusions, for us it remains to show that the 
infinitely many constraints on each Fj i implied by , viz. 


(sS’lsS’) 


VuG N, 


(41) 


where (. | .) now means inner product in £^(]R^"), do not impose impossible con- 

(n) 

sistency requirements. To show this, recall that the /g ' by definition satisfy 


[ ft~^^\vl,...,Vn+l)dVr^+l= (42) 

dK3 

which in view of implies that the hierarchy of the gj'^^ must satisfy 

3 

> Vn+l)dVn+l = ■■■,•«„) P 43 („+ i )_,. o , (43) 

i=l 

which is readily verified by explicit integration of 13511 . Thus, the constraints gu 
are automatically consistent, and this vindicates our initial assumption. 



5 Propagation of Chaos 

Setting n = 1 in (P|). and changing the time scale by setting r = |Eot, we recover 
113, with in place of /. However, (PI (or (0 for that matter) cannot be said 
to have been shown to be a kinetic equation yet. Note that propagation of chaos 
has not entered the derivation of (123. In fact, (123 for n = 1 , 2 ,... constitutes 
a “Fokker-Planck hierarchy” analogous to the the well-known Boltzmann, Landau 
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and Vlasov hierarchies which arise in the validation of kinetic theory |Spo91[RlIP94| 
using ensembles. In our case the hierarchy has the very simplifying feature that the 
n-th equation in the hierarchy is decoupled from the equation for the n + 1-th 
marginal. Since all the hierarchies used in the validation of kinetic theory are by 
construction linear^^ in the “vector” of the whenever one has a decoupling 

hierarchy one obtains closed linear equations for the . In particular, our equation 
with n = 1 is already a closed linear equation for However, at this point, 
any is still in general an ensemble superposition of states; in particular, 
still describes a statistical ensemble of pure states / with same mass, momentum, 
and energy. By ignoring this fact one can mislead oneself into thinking that 
with n = 1 and in place of / is already the kinetic equation we sought. 

The final step in extracting m as kinetic equation for the pure states involves 
the Hewitt-Savage |HeSa55| decomposition theorem. This theorem says that in the 
continuum limit any is a unique convex linear superposition of extremal (i.e. 
pure) n particle states, and that these pure states are products of n identical one- 
particle functions / evaluated at n generally different velocities. Each of the / in 
the support of the superposition measure represents the velocity density function of 
an actual individual member of the infinite statistical ensemble of infinitely-many- 
particles systems. In formulas, at r = 0 the initial data for /(”) read 

/(")(ui,...,u„;0) = (/o®"(ui,...,u„)), (44) 

where (.) is the Hewitt-Savage | HeSa55| ensemble decomposition measure on the 
space of initial velocity density functions /o of individual physical systems with same 
mass m.(/o)(= 1), momentum p{fo) — uq and energy e(/o) = cq = So + l'i^oP/2. To 
extend this representation to r > 0, let denote the one-parameter evolution 
semigroup for ((221), i.e. ... , t) = > '^n)- Noting now that 

the Hewitt-Savage measure is of course invariant under the evolution, and that by 
the linearity of (221 it commutes with the linear operator for all t ^ 0, it 

^^More precisely, they are only essentially linear, for the parameters eo and uq^ which also 
enter any of the other hierarchies whenever they describe ensembles of systems conserving mass, 
momentum, and energy, are all tied up with the initial conditions. 
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follows that at later times r > 0 the n point density of the ensemble is given by 


..., r) = v^)). (45) 

This so far simply states that, if the ensemble is initially a statistical mixture of pure 
states (product states), then at later times it is a statistical mixture of time-evolved 
initially pure states. Next we note that by inspection of 1221 it follows that 

C7('^)/o®”(^i,..., 7;„) = (C/W/o)®’^(t;i,...,^„), (46) 

viz. pure states evolve into pure states. Every factor f{vk',T) = C/r^Vo(^fc) solves 
m with r = obeying the desired conservation laws. At last one can legit¬ 
imately say that m has been derived as a full-fledged kinetic equation valid for 
almost every (w.r.t. (.)) individual member of the limiting ensemble. 

6 Summary and Outlook 

In summary, the diffusion equation on can be interpreted as the simplest 

“master equation” for an underlying iV-body Markov process with single-particle 
and pair terms. The N ^ oo limit for the marginal densities of solutions to the dif¬ 
fusion equation is well-defined and can be carried out explicitly. After invoking the 
Hewitt-Savage decomposition, the limit iV —> oo is seen to produce solutions of the 
“kinetic Fokker-Planck equation” describing individual isolated systems conserving 
mass, momentum, and energy. The Fokker-Planck equation is exactly solvable 
and displays correctly the qualitative behavior of a typical kinetic equation. In this 
sense, really can be regarded as the simplest example of a kinetic equation of 
the “diffusive” type, in the same family as, for instance, the much more complex 
Landau and Balescu-Lenard-Guernsey equations. 

Our work raises many new questions. 1) In particular, in Appendix Ab we have 
only written down the generator for the adjoint process of the underlying A^-particle 
Markov process; hence, what is the explicit characterization of this process? 2) A 
derivation of a kinetic equation a la Kac is an intermediate step towards a full 
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validation from some deterministic (Hamiltonian) microscopic model, which is in 
general a very difficult program, see the rigorous derivations of kinetic equations in 
|Spo91| [(11^941 . The substitute Markov process is usually chosen to preserve some 
of the essential features of the deterministic dynamics which (formally) leads to 
the same kinetic equation. Here we have only identified a stochastic model which 
leads to ©■ Villani’s work |Vil98| suggests that a deterministic model may exist 
which in the kinetic regime leads to 0. Can one indentify this model? 3) In this 
paper, we conveniently assumed that the Fourier coefficients ensure convergence of 
the marginal density functions in n and subsequently upgraded the regularity 
to Schwartz functions. What are the explicit conditions on the Fourier coefficients 
of the initial functions on which ensure convergence in n in Schwartz 

space, in some topology for measures? 4) Since the PDFs in our finite-iV Fokker- 
Planck hierarchy are already self-contained for each n (viz., they do not involve the 
usual coupling to the finite-?V corrections to the limiting evolutions can 

be studied in great detail; hence, for instance, how do the explicit corrections to 
propagation of chaos look? 5) We already mentioned in a footnote that the kinetic 
Fokker-Planck equation can easily be generalized to situations where the system 
is exposed to some external driving force by adding a forcing term. Can one de¬ 
rive this equation from some suitable ensemble of driven systems? Under which 
conditions do there exist stationary non-equilibrium states, and what are their sta¬ 
bility properties? 6) Finally, our derivation is only valid for the space-homogeneous 
Fokker-Planck equation without driving force term; hence, can one extend our 
derivation to obtain the space-inhomogeneous generalization of the kinetic Fokker- 
Planck equation, first without and then with driving force term? These are many 
interesting questions which should be answered in future works. 
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Appendix 

A. Two useful representations of the Laplacian on spheres 

Aa. Extrinsic representation in divergence form 

For the purpose of obtaining equations for the marginals by integrating 113, it 
is advantageous to express the Laplacian on the right-hand side in terms of the 
projection operator P^sN-i from to the fibers of the tangent bundle of the 
embedded manifold It is easy to verify |KiLa04| that 


= V • [Pjv 

“0 1^0 


r3N-4 

L-u.o,eo 




(47) 


In order to have an explicit expression for P„3iv-4 we introduce an orthogonal basis 
for the orthogonal complement of the tangent space to at T G C 


. Clearly, such orthogonal complement is spanned by the four vectors V and 
= (e^,..., e^), cr = 1,2, 3, where the are the standard unit vectors in 
The vectors are orthogonal to each other but not to V; projecting away the 
non-orthogonal component of V yields 

1 3 


(T—1 


V = V-U. 


(48) 


The vectors {V — U,Ei,E 2 ,Es} form the desired orthogonal basis; their magni¬ 
tudes are |Ecr| = VN and \ V — U\ = ■\/2Neo- Finally, 14711 becomes 


= dv ■ 

^'>‘^UQ,eQ 


IsAT — ^ (g) Ecr — 

\ cr=l 


2Neo 


{V-U)®iV-U)]dvF^^'> 


(49) 


Ab. Representation for the A-Body Markov Process 


In the main part of this paper we started from the diffusion equation on the man¬ 
ifold of A-body systems with same energy (per particle) eo and momen¬ 

tum (per particle) Mq; then took the limit N oo, obtaining the kinetic Fokker- 
Planck equation m, which rewrites into m in view of the conservation laws. The 
Laplace-Beltrami operator on is the generator of the adjoint semigroup of 
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the underlying stochastic Markov process that rules the microscopic dynamics of 
an individual A^-body system. Here we show that this generator can be written as 
a sum of single particle and two-particle operators, thus characterizing the Markov 
process as a mixture of individual stochastic motions and stochastic binary interac¬ 
tions. Moreover, we show that the binary particle operators are the only ones that 
do not vanish in the iV —> oo limit. This means that the kinetic Fokker-Planck 
equation can also be derived in terms of an A'^-body stochastic process with purely 
binary interactions, which is more satisfactory from a physical point of view. 

Recall that in section 2 we explained that can be identified with the 

sphere centered at the origin of (which itself is an affine linear sub¬ 
space of the space of all velocities, Recall that Ag.3 w-4 = As3w-4. Note 




the well-known representation 


_^ , V 2 

As3iv-4 = ^ [wkdwi - vuidwkj , (50) 

where Wk is the fc-th Cartesian component oi W € §3Af-4 ]^3(Af-i) 

in section 2 we used W G but note furthermore that the r.h.s. of is 

invariant under W ^ AW). Grouping the components of W into blocks of vectors 
Wk G A: = 1,..., N — 1, the r.h.s. of (15011 can be recast as 

N-l N-l 

A§3Ar-4 X! - {wk -dwi)^ 

fe=i 1=1 
l^k 

N-l 

-'^{wkX (51) 

fe=i 

containing one-body terms as well as binary terms. Note however that the first term 
in the binary sum is effectively a sum of two-body terms in disguise, which scale 
with factor N — 2 and thus survive in the limit N oo, while the true one-body sum 
(second line) drops out in that limit. This implies that the kinetic Fokker-Planck 
equation m can be derived from a master equation on which contains only 

the binary terms (first line) in (I51II . This in turn implies that m is the kinetic 
equation for an underlying system of N particles with stochastic pair interactions. 
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B. High-Dimension Asymptotics of Associated Legendre Func¬ 
tions 


In (TTHll the associated Legendre functions of degree s = 0,1,2,... and order r = 
0,..., s in g dimensions occur. They are defined on the interval [—1,1] and given by 


PIit;q) = Vq- 


■s+r ' 


9-1 


E(- 




(1 - t'^y+i 
4.J n{s-r-2iy.T{l + r + ^) 


(52) 


which differ from the PJ (t; q) in |Mul98| in their normalization. In our investigation, 


q = 3N — p and t = 


and we are interested in the limit N —>■ oo. 


The familiar asymptotics of Euler’s Gamma function gives us 


+ O . (53) 

L [a + X) V / 

for a; ^ 1. Applying this asymptotics with 2x = q — 1 = 3N — p— 1 and a = l + r to 
lO, we find that given p G N and w gM. (which implies N > max{p/3, (2eo)}), 
when ^ 1 we have 


pr 


W 




-.;3N-p 


1^1 

E (-1)' 

1=0 


n{s-r-2l)\ 




s—r—2l 


(54) 


By comparing with the formula for the Hermite polynomial of degree k on K, 

[IJ 

H,{x) = E(-l)'+'=yj^^^^(2a.)'=-^', (55) 

we see that, given p G N and w G M, we have 

+ °{ 7 n ) 

when A ^ 1. Hence, for all fixed p we now find that pointwise for any w G K., 

(ys- “ - ”) = (-^)'" (\/5“) 

where again it is understood that A > max{p/3, w'^/ (2eo)} in the expression under 
the limit in the left-hand side. Equation (I35II in the main text follows. 
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